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Abstract 

Axial  double-ball  load  tests  were  made  on  freshwater  ice  and  first-year  and 
multiyear  sea  ice.  From  this  simple  test  method  the  apparent  unconfined 
compressive  strength  of  the  ice  was  determined.  These  strength  results  were 
compared  with  those  obtained  from  the  complex  and  costly  uniaxial  unconfined 
compression  test  made  on  similar  ice  at  a  strain  rate  of  10-3  s-1 .  The  scatter 
in  the  test  data  and  the  average  ice  strength  obtained  from  both  test  methods 
were  similar.  The  findings  indicate  that  the  expedient  axial  double-ball  load  test 
is  well  suited  for  determining  the  unconfined  compressive  strength  of  ice, 
especially  in  the  field  where  the  demanding  sample  preparation  requirements 
needed  for  unconfined  compression  test  samples  cannot  be  met. 


For  conversion  of  Si  metric  units  to  U.S./British  customary  units  of  measurement 
consult  Standard  Practice  for  Use  of  the  International  System  of  Units  (SI),  ASTM 
Standard  E380-89a,  published  by  the  American  Society  for  Testing  and  Mater¬ 
ials,  1916  Race  St.,  Philadelphia,  Pa.  19103. 
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Axial  Double-Ball  Test  Versus  the 
Uniaxial  Unconfined  Compression  Test  for 
Measuring  the  Compressive  Strength  of 
Freshwater  and  Sea  Ice 

AUSTIN  KOVACS 


INTRODUCTION 

Ice  formed  in  nature  is  almost  always  structur¬ 
ally  anisotropic.  An  exception  is  the  equiaxed  gran¬ 
ular  ice  found  at  select  depths  in  polar  ice  caps, 
where  the  c-axes  of  the  ice  crystals  are  randomly 
oriented.  This  ice  has  been  found  to  exhibit  the 
same  strength  properties  irrespective  of  the  direc¬ 
tion  of  the  applied  load  (Kovacs  et  al.  1969).  Fresh¬ 
water  ice  is  often  considered  to  be  a  single-phase 
material;  however,  it  frequently  contains  numer¬ 
ous  gas  bubbles  of  varying  size  that  affect  ice 
strength.  Sea  ice  is  truly  a  multiphase  material 
composed  of  solid  freshwater  ice,  liquid  brine  and 
gas.  Like  freshwater  ice  the  sea  ice  structure  can 
consist  of  granular  and  columnar  crystals,  which 
can  vary  in  width,  length  and  structure  with  ice 
sheet  depth  .Most  natural  ice  is  heterogeneous  and 
anisotropic;  therefore,  the  structure  of  the  ice  and 
the  direction  of  the  load  applied  to  it,  as  well  as  its 
temperature  and  the  rate  of  load  application,  must 
be  factored  into  the  interpretation  of  load  test 
results. 

Many  test  techniques  have  been  applied  to  un¬ 
derstanding  the  strength  characteristics  of  ice,  with 
mixed  results.  The  most-performed  test  is  the  uniax¬ 
ial  unconfined  compression  test.  This  is  a  decep¬ 
tively  simple  test  that  has  not  always  been  done 
"correctly"  because  of  the  difficult  and  labor-in¬ 
tensive  sample  preparation  required.  In  addition, 
a  stiff  testing  machine  with  sample  load  or  strain 
rate  feedback  control  is  needed.  Another  issue  still 
being  reviewed  is  the  preferred  contact  that  should 
exist  between  the  test  sample  and  the  loading 


apparatus.  To  help  resolve  some  of  the  problems 
and  to  attempt  to  develop  a  standard  test  proce¬ 
dure,  the  International  Association  of  Hydraulic 
Research  Committee  on  Ice  Problems  devised 
guidelines  for  performing  the  uniaxial  unconfined 
compression  test  on  ice  (Schwarz  et  al.  1981).  Nev¬ 
ertheless,  sample  size  and  platen  contact  condi¬ 
tions  are  still  a  topic  of  review  (Kuehn  et  al.  1992), 
and  the  equipment  required  to  make  high-quality 
ice  samples  and  the  heavy  loading  apparatus  need¬ 
ed  to  perform  the  uniaxial  unconfined  compres¬ 
sion  test  pose  logistic  problems  when  field  testing 
is  desired. 

To  circumvent  the  demanding  sample  prepara¬ 
tion  requirements  and  heavy  equipment  needed 
to  perform  uniaxial  unconfined  compression  tests 
in  the  field,  Kovacs  (1978)  explored  the  feasibility 
of  using  a  simple  lightweight  axial  double-ball 
load  test  device  for  determining  the  unconfined 
compressive  strength  <jc  of  snow  and  ice  in  the 
Antarctic.  This  evaluation  led  to  another  study  in 
which  sea  ice  strength  was  evaluated  vs.  loading 
ball  size,  temperature  and  ice  density  (Kovacs 
1985).  From  this  study  the  unconfined  compres¬ 
sive  strength  determined  with  the  use  of  16-mm- 
diameter  loading  balls  was  found  to  agree  with  the 
uniaxial  unconfined  compressivestrength  obtained 
at  a  strain  rate  of  10-3  s-1. 


AXIAL  DOUBLE-BALL  TEST 

The  CRREL  axial  double-ball  (axial  DB)  test 
system  consists  of  a  test  frame,  a  manually  operat- 
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Figure  1.  Axial  DB  test  apparatus  in  use  at  CRREL. 


ed  hydraulic  pump  and  ram,  and  a  digital  load 
indicator  (Fig.  1).  The  test  frame  consists  of  two 
upright  stainless  steel  plates.  The  distance  be¬ 
tween  the  plates  is  adjustable  to  accommodate  ice 
samples  of  different  lengths.  After  adjustment  the 
plates  are  locked  parallel  into  a  rigid  structure  by 
tightening  bolts.  To  the  inside  of  one  plate  is  mount¬ 
ed  a  hydraulic  ram,  and  on  the  opposite  plate  is 


fastened  a  load  cell.  On  the  load  cell  and  on  the 
hydraulic  ram  are  mounted  load  balls.  To  support 
the  ice  sample,  two  parallel  rods,  spaced  about  3 
cm  apart,  extend  between  the  plates.  These  rods 
can  be  adjusted  up  or  down  to  accommodate  ice 
samples  of  different  diameters  to  ensure  that  the 
load  balls  are  in  line  with  the  axis  of  the  ice  sample. 
The  hydraulic  hand  pump  w  as  sized  to  ensure  that 
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Figure  2.  Axial  DB  test  jig  with  an  LVDT 
configured  to  measure  hydraulic  ram  travel 
vs.  time.  The  arrow  points  to  the  LVDT 
motion  shaft.  The  test  sample  is  sea  ice. 


Figure  3.  Logging  sea  ice  core  in  the  field. 
The  core  length  is  first  measured.  Small  holes 
are  then  drilled  about  2-3  cm  deep  into  the 
core  side  at  20-cm  increments,  and  a  ther¬ 
mistor  is  inserted  in  each  hole  to  determine 
the  in  situ  temperature  of  the  ice  sheet.  The 
core  is  next  cut  into  lengths  using  the  cutoff 
saw  on  the  right,  and  the  axial  length  of  the 
ice  cylinder  is  determined  (in  this  case  using 
a  dial  gauge  measurement  jig  located  beside 
the  notebook).  The  sample  is  nextweighedon 
an  electronic  balance  and  tested  in  the  axial 
DB  system.  The fractured  ice  pieces  are  then 
placed  in  plastic  bags  for  later  melting  and 
salinity  measurements. 


one  swift  downward  stroke  of  the  pump  handle 
would  quickly  drive  the  ball,  mounted  on  the  face 
of  the  hydraulic  ram,  sufficiently  deep  into  the  ice 
to  cause  it  to  split.  For  several  tests  a  linear  poten¬ 
tiometer  was  installed  for  monitoring  the  ram 
travel  and  the  time  to  peak  load  (Fig.  2). 

Test  sample  preparation  is  straightforward. 
Previous  testing  revealed  that  for  ice  the  sample 
length  must  be  1.05-1.1  times  the  sample  diame¬ 
ter.  The  former  length  is  preferred.  Ice  samples 
should  have  reasonably  parallel  ends,  which  need 
not  be  smooth.  This  may  be  achieved  with  use  of  a 
miter  box  or,  as  we  prefer,  with  the  use  of  a  carpen¬ 
ter's  power  cutoff  saw  (Fig.  3) .  With  this  device,  ice 
core  may  be  quickly  cut  to  length  in  the  field  and 
then  weighed  on  a  small  battery-powered  elec¬ 
tronic  scale.  The  sample  diameter  once  determined, 
using  a  pi  tape,  should  remain  constant  for  the  core 


barrel  used.  The  sample  length  may  be  deter¬ 
mined  with  a  caliper.  However,  we  prefer  a  digital 
electronic  displacement  measurement  device  for 
measuring  the  sample  length  at  the  axis  of  the  ice 
core.  The  simple  measurement  jig  made  for  field 
use  is  shown  in  Figure  4.  This  measurement  tech¬ 
nique  is  very  accurate  and  compensates  for  sam¬ 
ple  ends  not  being  parallel  should  that  occur  when 
a  hand  saw  is  used  to  cut  core  to  length.  An  accu¬ 
rate  length  measurement  is  most  needed  for  the 
determination  of  sample  volume  but  also  for  use 
in  the  axial  DB  strength  equation. 

The  test  procedure  is  also  straightforward.  An 
ice  sample  is  set  onto  the  support  rods,  the  ram  is 
moved  forward  until  the  balls  contact  the  ice,  and 
the  ram  is  then  propelled  forward  with  one  down- 
stroke  of  the  hydraulic  pump  handle.  When  the 
balls  are  forced  against  the  ice,  stresses  develop  in 
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b.  Measuring  the  core  diameter. 

Figure4.  Digital  measuring  device.  The  cable  provides  power  to  a  small 
blanket  heater  needed  to  prevent  the  crystal  display  from  freezing. 


Load  Ball 


Sample  deaves 

when  crack  spontaneously  progresses 
through  critically  stressed  region. 
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Figure  5.  Failure  concept  for  ball  loading  of  a 
test  sample  as  postulated  by  Reichmuth  (1968). 


the  ice  at  the  contact  zone.  In  the  freshwater  ice 
sample  shown  in  Figure  1,  cracks  have  developed 
in  the  stress  zone,  giving  the  ice  a  cloudy  white 
appearance  at  each  end.  The  crack  formation  se- 
quer  ~e  has  been  described  by  Reichmuth  (1968). 
He  speculated  that  en  echelon  cracks  first  form  in 
the  stressed  zone  (Fig.  5).  As  the  cracks  multiply, 
they  coalesce  and  form  sliplines.  Sample  failure 


J- 


occurs  when  a  critical  stress  level  is  reached,  which 
most  often  causes  a  crack  to  propagate  along  the 
axis  of  the  sample.  The  peak  force  at  failure,  as 
sensed  by  the  load  cell,  is  displayed  on  the  high¬ 
speed  digital  peak-load  indicator.  The  unique  fea¬ 
ture  of  this  test  procedure  is  that  ice  may  be  tested 
shortly  after  being  removed  from  the  core  barrel. 
Failed  ice  samples  typically  split  in  two  (Fig.  6)  or 
three  pieces  (Fig.  7),  but  four  pieces  sometimes 
occur.  Spalling  also  occurs  in  conjunction  with 
axial  splitting  and  at  times  as  the  sole  fracture  (Fig. 
8).  The  latter  event  renders  the  test  unacceptable. 


a.  Sea  ice.  b.  Freshwater  ice. 

Figure  6.  Example  of  failed  ice  samples  loaded  in  the  axial  DB  test  system.  Note  the  interesting  nonlinear  split  in 
the  sea  ice  sample  and  the  lacy  crack  “pipe"  visible  at  the  center  of  the  right  freshwater  ice  piece.  "Cloud"  cracks  are 
also  visible  in  the  freshwater  ice. 


Figure  7.  Multiyear  sea  ice  test  sample  that  Figure  8.  Multiyear  sea  ice  sample  that  failed  by  spalling, 
split  into  three  segments  at  failure. 
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Figure  9.  Shape  factor  vs.  ice  core  diameter.  (After  Kovacs  1978.) 


This  type  of  failure  may  be  the  result  of  a  slip-line 
crack,  shown  in  Figure  5,  propagating  laterally. 

The  equation  for  determining  oc  from  the  axial 
DB  test  is  (Kovacs  1978) 


where  K-  a  shape  constant 
P  =  peak  force  (lbf) 

L  =  sample  length  (in.). 

The  shape  constant  varies  with  sample  diame¬ 
ter  and  is  determined  as  shown  in  Figure  9.  K  is 
used  to  normalize  the  peak  force  resulting  from 
the  testing  of  samples  with  different  diameters. 

ICE  TESTED 

Laboratory-grown  freshwater  ice  and  natural 
first-year  and  multiyear  sea  ice  were  tested.  Labo¬ 
ratory  ice  was  made  at  CRREL  using  a  procedure 
described  by  Cole  (1979)  and  at  the  Ice  Research 
Laboratory  (IRL)  of  the  Thayer  School  of  Engineer¬ 
ing,  Dartmouth  College,  using  a  modified  Cole 
procedure  (Schulson  1990).  The  CRREL-grown 
ice,  made  in  ~70-mm-diameter  molds,  had  a  mean 
density  of  0.904  Mg/m3  and  a  related  porosity  of 
-1.4%.  The  equiaxed  grains  had  a  mean  diameter 
of  3.5  mm.  The  same  ice  was  being  used  at  IRL  in 
an  assessment  of  the  effect  of  grain  size  on  the 
uniaxial  oc  (Cannon  1985,  Schulson  1990).  This  ice 
was  made  in  -82-mm-diameter  molds  and  con¬ 
tained  grains  having  a  mean  size  of  5  mm.  The  ice 
had  a  preferred  higher  density  of  -0.912  Mg/m3 
and  a  lower  porosity  of  -0.6%. 


The  multiyear  sea  ice  tested  with  the  axial  DB 
loading  system  was  collected  in  the  Beaufort  Sea 
by  Cox  et  al.  (1984).  They  made  uniaxial  uncon¬ 
fined  compression  tests  on  the  same  source  ice  at 
CRREL  using  right  cylinders  having  a  machined 
diameter  of  102  mm.  We  tested  samples  at  CRREL 
using  the  as-cored  diameter  of  -105  mm. 

The  first-year  sea  ice  cores  were  obtained  from 
large  1 .2-m-square  blocks  of  ice  removed  from  the 
1.8-m-thick  sea  ice  in  Stefanson  Sound,  Alaska. 
The  ice  cores  were  mined  from  the  ice  block  at  the 
same  depth  in  both  the  vertical  (0°)  and  horizontal 
(90°)  plane  of  the  ice  sheet.  In  the  latter  the  ice  cores 
were  taken  parallel  90°-0°  and  perpendicular  90°- 
90°  to  the  well-aligned  horizontal  c-axes  of  the 
columnar  ice  crystals.  The  mean  crystal  size  was 
about  11  mm.  The  77-mm-diameter  ice  cores  were 
cut  to  length  and  tested  in  a  shed  on  shore. 

Vertical  ice  cores  were  also  obtained  from  the 
first-year  sea  ice  in  a  refrozen  melt  pool  in  a  mul¬ 
tiyear  floe  and  from  the  fast  ice  in  Stefanson  Sound. 

TEST  RESULTS 
Freshwater  ice 

Six  ice  cylinders  were  provided  from  the  IRL  at 
Dartmouth  College.  From  these  cylinders,  12  sam¬ 
ples  were  cut.  Each  had  a  length  of  -102  mm  and 
a  diameter  of  91.4  mm.  One  of  the  tested  samples 
is  shown  in  Figure  10.  Note  the  crushed  zone  at 
each  end  of  the  ice  cylinder  and  what  appears  to  be 
eti  echelon  cracks  extending  below  the  top  crushed 
zone.  The  fractured  ice  surface  is  laced  with  crack 
damage,  and  there  is  a  vertical  "pipe"  of  highly 
fractured  ice  containing  a  fine  lacework  of  cracks. 
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Figure  l 0.  Example  of  failed freshwater  ice  show¬ 
ing  an  axial  "pipe"  of  concentrated  fine  cracks, 
en  echelon  cracks  (onein front  of  the  arrow)  and 
numerous  "cloud"  cracks  that  may  have  formed 
at  crystal  boundaries  within  the  ice. 

Table  1.  Test  data  for  5-mm  grain  size  freshwater  ice  at  -8°C. 


Young's  Shear  Volume 


Sample 

no. 

Specific 

gravity 

Velocity  (m/s) 

P  wave  S  wave 

Poisson's 

ratio 

modulus 
( CPa ) 

modulus 

(GPa) 

rigidity 

(GPa) 

Oc 

(MPa) 

1L 

0.910 

3845 

1728 

0.374 

7.5 

2.7 

9.8 

7.81 

1R 

0.912 

3817 

1629 

0.389 

6.7 

2.4 

10.0 

7.03 

2L 

0.913 

3799 

1697 

0.375 

7.2 

2.6 

9.7 

6.44 

2R 

0.912 

3822 

2001 

0.311 

9.6 

3.6 

8.5 

4.15 

3L 

0.909 

3811 

1704 

0.375 

7.3 

2.6 

9.7 

4.66 

3R 

0.913 

3784 

1680 

0.377 

7.1 

2.6 

9.6 

5.18 

4L 

0.915 

3761 

1689 

0.374 

7.2 

2.6 

9.5 

5.00 

4R 

0.914 

3781 

1660 

0.381 

7.0 

2.8 

10.8 

5.12 

5L 

0.908 

3812 

1703 

0.375 

7.2 

2.6 

9.7 

4.42 

5R 

0.913 

3818 

1854 

0.346 

8.4 

3.1 

9.1 

8.06 

6L 

0.914 

3819 

1984 

0.315 

9.5 

3.6 

8.5 

5.54 

6R 

0.913 

3812 

1993 

0.284 

9.5 

4.6 

8.5 

5.29 

Avg 

0.912 

3807 

1765 

0.358 

7.9 

2.9 

9.5 

5.72 

The  -8°C  test  results  are  listed  in  Table  1.  The 
average  oc  strength  was  5.72  MPa,  which  needs  to 
be  adjusted  to  -10°C,  the  temperature  at  which  10 
right  cylinders  of  the  same  ice  type  were  tested,  at 
the  IRL,  in  uniaxial  unconfined  compression.* 
Schulson  stated  that  the  IRL  tests  gave  a  mean  ice 
strength  of  5.8  MPa  at  a  strain  rate  of  10-3  s-1. 

To  adjust  oc  for  temperature,  a  correction  factor 
based  on  the  analysis  of  Kovacs  et  al.  (1977)  was 
used.  They  found  that  oc  for  snow  and  ice  changed 
at  a  rate  of  -0.075  MPa/°C.  This  value  is  in  good 
agreement  with  the  work  of  Butkovich  (1954)  on 
lake  ice,  Wolfe  and  Thiem  (1964)  on  river  ice, 
Kovacs  (1978)  on  multiyear  sea  ice  and  Brown  and 


*  Personal  communication.  Dr.  Erland  Schulson,  Dartmouth 
College. 


McKittrich  (1992)  on  a  single  crystal  of  freshwater 
ice.  These  authors  found  oc  to  increase,  at  temper¬ 
atures  below  about  -10°C,  at  a  rate  of  0.074, 0.08, 
0.07  and  0.065 MPa/°C,  respectively.  Bender  (1957) 
proposed  the  following  empirical  expression  to 
account  for  the  change  in  oc  with  temperature: 

oc2/acl=m/\nr6  ro 

where  Cc,  and  CTq  are  the  oc  strengths  at  tempera¬ 
tures  Ti  and  Tz,  respectively.  Based  on  the  findings 
of  the  above-referenced  authors,  this  equation  pro¬ 
vides  a  strength  correction  that  is  too  high.  Chang¬ 
ing  the  exponent  to  0.13  would  bring  the  Bender 
equation  into  good  agreement  with  a  oc  strength 
change  vs.  temperature  of  0.075  MPa/°C. 

When  this  temperature  correction  is  applied  to 
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the  average  axial  DB  oc  value  of  5.72  MPa  at  -8°C, 
it  increases  to  5.87  MPa  at  -10°C.  Clearly  the  aver¬ 
age  axial  DB  stiength  is  in  excellent  agreement 
with  the  average  uniaxial  oc  value  (5.8  MPa)  pro¬ 
vided  by  Schulson  for  the  1RL  tests.  The  two  inde¬ 
pendent  test  techniques  on  the  same  ice  type  show 
that  the  simple  axial  DB  test  provides  comparable 
oc  values  to  those  obtained  in  the  laboratory  at 
lfr^s-1  using  highly  sophisticated  testing  equip¬ 
ment  and  sample  preparation  techniques. 

Prior  to  axial  DB  testing,  a  sonic  velocity  meter 
was  used  to  measure  the  100- kHz  compression 
wave  velocity  Vp  and  shear  wave  velocity  Vs  in 
each  sample.  The  objective  of  making  the  Vp  and 
Vs  wave  determinations  was  to  see  if  a  correlation 
existed  between  oc  and  the  velocity  values.  With 
the  Vs  and  Vp  measurements  and  the  previously 
determined  ice  density  (specific  gravity  y),  the 


parameters  shown  in  the  following  flow  chart 
were  calculated: 


where  v  =  Poisson's  ratio 
S  =  shear  modulus 
£  =  Young's  modulus 
R  =  volume  rigidity. 

The  average  dynamic  Young's  and  shear  mod¬ 
uli  given  in  Table  1  for  the  axial  DB  tested  ice  were 


a.  Young's  modulus  vs.  density. 


b.  Shear  modulus  vs.  density. 


Figure  11.  Dynamic  modulus  test  results  for 
polar  fim  and  ice  (Smith  1965)  and  the 
average  test  data  from  this  study. 
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Table  2.  Axial  double-point  unconfined  compressive  strength  vs.  density  of  freshwater  ice 
with  a  grain  size  of  3.5  mm. 


-W°C  tests 

-20*C  tests 

-10‘ C  tests 

-20t  tests 

Sample 

no. 

Density 

(Mg/m3) 

Oc 

(MPa) 

Sample 

no. 

Density 

(Mg/m3) 

oc 

(MPa) 

Sample 

no. 

Density 

(Mg/m3) 

Oc 

(MPa) 

Sample 

no. 

Density 

(Mg/m3) 

Oc 

(MPa) 

1A 

0.909 

7.40 

IB 

0.912 

10.78 

11A 

0.906 

9.23 

11B 

0.905 

7.51 

2A 

0.907 

7.71 

2B 

0.907 

7.69 

12A 

0.905 

8.47 

12B 

0.904 

8.79 

3A 

0.891 

7.57 

3B 

0.890 

8.37 

13A 

0.904 

7.23 

13B 

0.907 

9.19 

4A 

0.906 

7.47 

4B 

0.908 

8.46 

14A 

0.893 

8.05 

14B 

0.875 

7.45 

5A 

0.905 

7.82 

5B 

0.907 

7.84 

15A 

0.908 

6.99 

15B 

0.909 

9.02 

6A 

0.902 

7.72 

6B 

0.909 

7.20 

16A 

0.904 

7.64 

16B 

0.905 

7.36 

7A 

0.903 

7.41 

7B 

0.907 

7.84 

17A 

0.908 

8.63 

17B 

0.908 

10.47 

8A 

0.904 

7.43 

8B 

0.909 

8.92 

18A 

0.896 

8.00 

18B 

0.898 

9.48 

9A 

0.907 

6.95 

9B 

0.906 

8.97 

19A 

0.886 

7.73 

19B 

0.887 

7.42 

10A 

0.898 

7.99 

10B 

0.897 

7.58 

Avr 

0.902 

7.76 

Avg 

0.903 

8.44 

i 


Figure  12.  Freshwater  ice  unconfined 
compressive  strength  vs.  grain  size  at 
-10°C  and  a  strain  rate  of  lO'3  s1. 


7.9  and  2.9  GPa,  respectively.  These  average  values 
compare  extremely  well  with  the  dynamic  modu¬ 
lus  test  results  of  Smith  (1965),  who  tested  fim  and 
ice  at  Camp  Century,  Greenland,  between  -12°  and 
-15°C  (Fig.  11).  Also,  the  average  value  for  Pois¬ 
son's  ratio  in  Table  1  agrees  with  the  findings  of 
Smith.  No  temperature  correction  was  applied  to 
the  modulus  data,  as  Mellor  (1983)  has  emphasized 
that  the  modulus  of  ice  is  not  sensitive  to  tempera¬ 
ture  in  the  range  of  these  tests.  No  correlation  was 
found  between  oc  and  the  Vp  and  Vs  wave  determi¬ 
nations  made  on  the  ice,  which  had  a  specific  grav¬ 
ity  of  -0.91. 

The  axial  DB  test  results  made  at  -10°  and  -20°C 
on  the  3.5-mm-grain-size  freshwater  ice  made  at 
CRREL  are  listed  in  Table  2.  The  average  oc  values 
at  -10°  and  -20°C  were  7.76  and  8.44  MPa,  respec¬ 
tively.  The  change  in  ice  strength  from  -10°  to 
-20°C  is  0.068  MPa/°C,  which  is  in  agreement  with 
the  ctc  temperature  correction  previously  discussed. 

To  assess  how  well  the  temperature-corrected 
axial  DB  oc  values  agree  with  uniaxial  unconfined 


compressive  strength  tests  made  at  10“3  s-1,  the 
axial  DB  oc  data  are  compared  with  -10°C  uniaxial 
ac  data  for  ice  of  the  same  grain  size.  In  Figure  12 
the  uniaxial  oc  test  results  of  the  authors  listed  and 
the  average  axial  DB  oc  values  obtained  for  the  ice 
with  3.5-  and  5-mm  grain  sizes  are  shown.  The 
axial  DB  values  fall  nicely  beside  the  appropriately 
weighted  regression  curve  passing  through  the 
data.  The  outlier  at  the  far  right  was  not  included 
in  the  regression  analysis. 

The  data  point  shown  as  Camp  Century  in 
Figure  12  was  determined  as  follows.  In  1966  Gow 
(1975)  determined  the  grain  size  of  the  fim  and  ice 
in  the  inclined  drift  at  Camp  Century,  Greenland. 
This  sloping  passageway  (Fig.  13)  extended  to  a 
depth  of  100  m  below  the  surface,  where  the  ice 
density  reached  a  density  of  0.89  Mg/m3  at  an 
ambient  temperature  of  -24°C.  We  obtained  his 
data*  and  analyzed  it  as  shown  in  Figure  14.  When 
the  regression  curve  in  Figure  14  is  extrapolated  to 


*  Personal  communication.  Dr.  Anthony  Gow,  CRREL. 
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,  Unconfined  Compressive  Strength  (MPa)  Q> Grain  S“®  (mm) 


Figure  13. 1966  view  looking  down  the 
inclined  drift  at  Camp  Century ,  Green¬ 
land,  where  homogenous  ice  samples  were 
obtained.  The  pipe  and  metal  objects  beside 
the  pipe  were  part  of  two  conveyor  systems 
used  to  remove  material  during  mining  of 
the  drift. 


p.  Density  (Mg/m3)  Gow,  private  communication.) 


Figure  15.  Camp  Century,  Greenland,  fim 
and  ice  unconfined  compression  strength 
os.  density  at  -25  °C  and  a  strain  rate  of 
~20"3  s'1.  (From  Kovacs  et  al.  1969.) 
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0.92  Mg/m3,  the  density  of  pure  ice  at  -24°C,  we 
obtain  a  grain  size  of  1.76  mm. 

In  1967,  uniaxial  unconfined  compression  tests 
were  made  in  the  Camp  Century  inclined  drift 
(Kovacs  et  al.  1969).  These  tests  were  made  at  a 
strain  rate  of  about  10-3  s_1  and  a  temperature  of 
-25°C.  These  test  data  are  replotted  in  Figure  15. 
The  regression  curve  passing  through  the  test  re¬ 
sults,  when  extrapolated  to  a  density  of  0.92  Mg/ 
m3,  gives  a  uniaxial  unconfined  compressive 
strength  of  11.11  MPa.  Using  the  oc  temperature 
correction  of  0.075  MPa/°C,  this  uniaxial  oc  value 
reduces  to  9.99  MPa  at  -10°C,  where  the  ice  on 
warming  would  now  have  a  density  of  -0.917 
Mg/m3.  This  uniaxial  oc  value  is  plotted  vs.  the 
above-determined  grain  size  of  1 .76  mm  in  Figure 
12. 

There  is  considerable  data  scatter  in  Figure  15. 
However,  such  scatter  is  common  with  uniaxial 
unconfined  compression  tests.  There  are  numer¬ 
ous  reasons  for  this  scatter,  including  sample  prep¬ 
aration  anomalies,  internal  flaws,  type  of  ice- 
machine  contact,  variations  in  sample  alignment 
within  the  testing  machine  as  well  as  ice  density 
(i.e.,  porosity)  variation.  Gow  (1975)  pointed  out 
that  bubbles  represent  defects  in  the  ice  that  "could 
be  expected  to  influence  the  strength  of  the  ice 
until  they  are  entirely  eliminated"  as  they  act  as 
stress  concentrators  that  could  lead  to  a  lower 
peak  strength.  In  any  event,  what  Figure  15  does 
indicate  is  the  need  to  run  multiple  tests  on  the 
same  ice  in  order  to  obtain  a  representative  mean 
peak  strength  value.  In  our  previous  work  we 
attempted  to  use  10  or  more  samples  to  obtain  an 
average  value  (Kovacs  et  al.  1969,  1977,  Kovacs 
1985).  This  minimum  number  is  in  agreement  with 
the  work  of  Yamaguchi  (1970),  who  determined 


that  a  minimum  of  10  samples  are  needed,  and 
more  recently  by  Kirk  (1989),  who  indicated  that 
12  samples  is  a  reasonable  test  number.  This  many 
samples  can  be  difficult  to  obtain  and  transport  to 
the  testing  facility.  The  ice  at  Camp  Century  was 
unique  in  that  at  a  given  depth  one  could  extract 
endless  samples  of  the  same  density  and  grain 
structure.  Unfortunately  this  resource  is  no  longer 
accessible.  In  principle,  laboratory-grown  ice  can 
also  be  replicated  to  a  high  degree  of  structural  and 
density  conformity.  However,  first-year  sea  ice 
does  not  lend  itself  to  repetitive  sample  conformi¬ 
ty  due  to  a  wide  variation  in  this  material's  compo¬ 
sition  and  structure  within  the  ice  sheet  (Weeks 
and  Ackley  1982). 

Multiyear  sea  ice 

The  physical  and  mechanical  properties  of  the 
multiyear  sea  ice  tested  are  listed  in  Table  3.  Also 
shown  is  the  time  to  failure  and  the  total  distance 
that  the  two  load  balls  penetrated  the  ice  at  peak 
loading.  The  time  to  failure  in  the  axial  DB  test  is 
more  than  an  order  of  magnitude  shorter  than  the 
typical  times  to  failure  in  a  standard  uniaxial  un¬ 
confined  compression  test  run  at  a  10~3  s_1  strain 
rate. 

The  axial  DB  oc  values  listed  in  Table  3  show 
considerable  scatter.  This  is  common  with  strength 
data  obtained  from  uniaxial  unconfined  compres¬ 
sion  testing  of  multiyear  sea  ice.  To  illustrate  this 
and  compare  our  results  with  the  uniaxial  oc  val¬ 
ues  obtained  by  Cox  et  al.  (1984)  at  -10°C  on  a 
universal  electrohydraulic  testing  machine,  we 
first  applied  the  temperature  correction  of  0.075 
MPa/°C  to  the  axial  DB  oc  data  obtained  at  -5°C. 
Both  data  sets  are  shown  in  Figure  16.  The  axial  DB 
oc  test  results  fall  well  within  those  obtained  by 


Table  3.  Test  data  for  multiyear  sea  ice  at  -5°C. 


Sample 

no. 

Brine 

vol. 

(%*) 

Air 

vol. 

(%o) 

Porosity 

(%o) 

Bulk 

density 

(Mg/m3) 

Ice 

density 

(Mg/m3) 

Failure 

load 

(kg) 

Time 
to  failure 
(ms) 

Ball 

travel 

(mm) 

Oc 

(MPa) 

RIA-a 

32.5 

9.4 

41.9 

0.913 

0.888 

254 

96 

7.1 

7.21 

RIA-b 

32.5 

9.4 

41.9 

0.908 

0.887 

211 

91 

6.0 

5.09 

RID 

4.8 

21.2 

26.0 

0.899 

0.894 

182 

80 

2.5 

6.05 

C8-a 

15.7 

8.8 

24.5 

0.912 

0.895 

200 

50 

2.8 

5.59 

C8-b 

15.7 

7.7 

23.4 

0.913 

0.896 

225 

35 

3.4 

6.19 

C9-a 

9.7 

15.4 

25.2 

0.905 

0.895 

202 

40 

6.9 

5.74 

C9-b 

19.6 

13.8 

33.4 

0.908 

0.887 

209 

90 

5.9 

5.77 

C9< 

16.7 

6.8 

23.5 

0.914 

0.896 

216 

40 

5.4 

5.96 

Cll-a 

9.8 

12.2 

22.0 

0.908 

0.898 

211 

72 

3.7 

6.00 

Cll-b 

12.7 

10.5 

23.2 

0.910 

0.8% 

257 

42 

3.2 

7.11 

Cll-c 

12.7 

10.5 

23.2 

0.910 

0.8% 

209 

25 

4.1 

5.89 

Cll-d 

17.7 

8.0 

25.7 

0.913 

0.894 

254 

110 

— 

6.99 

Cll-e 

17.7 

9.1 

26.8 

0.912 

0.893 

204 

52 

3.6 

5.11 

Avg 

16.8 

11.0 

27.8 

0.910 

0.893 

218 

63 

4.6 

6.05 
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Figure  16.  Multiyear  sea  ice  unconfined  com¬ 
pressive  strength  vs.  porosity  at  -10°C  and  a 
strain  rate  of  10~3  s_1. 


Cox  et  al.  (1984).  However,  unlike  their  tests  the 
axial  DB  test  does  not  give  stress-strain  data  from 
which  an  estimate  of  Young's  modulus  or  Pois¬ 
son's  ratio  can  be  made. 

The  regression  curve  passing  through  the  data 
in  Figure  16  indicates  that  the  strength  of  multi¬ 
year  sea  ice,  as  with  all  materials,  is  governed  by  its 
porosity.  However,  many  more  test  results  are 
required  to  better  define  this  trend,  not  only  at  this 
strain  rate  but  at  other  strain  rates  and  tempera¬ 
tures  as  well. 

First-year  sea  ice 

The  strength  of  first-year  sea  ice  is  more  diffi¬ 
cult  to  categorize  because  of  its  wide  structure 
anisotropy  and  inhomogeneity.  Brine  volume,  ice 
density  and  porosity  are  the  most-used  sea  ice 
physical  properties  against  which  sea  ice  strength 
hasbeen  correlated .  "Hie  latter  is  preferred  because 
porosity  controls  the  strength  of  a  material  when 
all  other  conditions  are  equal. 

Vertical  and  horizontal  axial  DB  load  tests  were 
made  on  the  sea  ice  from  Steffanson  Sound.  As 
previously  mentioned,  this  ice  contained  ice  crys¬ 
tals  about  11  mm  in  width.  The  c-axes  were  highly 
aligned  in  the  horizontal  plane.  Tests  were  made  at 
-5°  and  -16°C,  as  shown  in  Tables  4-6.  Little  has 
been  published  on  sea  ice  strength  vs.  temperature 
and  what  has  is  conflicting  (e.g.  Schwarz  [1971] 
showed  a  high  oc  change  with  temperature,  while 
Brown  and  McKittrich  [1992]  show  the  opposite). 
Therefore,  the  ac  data  are  referenced  to  the  ice 
porosity,  a  common  procedure  in  sea  ice  strength 
studies. 

Axial  DB  load  tests  were  applied  perpendicular 
to  the  predominant  c-axis  orientation,  or  in  the 


convention  of  Peyton  (1966),  at  90°-90°,  and  paral¬ 
lel  to  the  preferred  c-axis  alignment  at  90°-0°.  The 
direction  of  the  applied  load  was  a  visual  estimate. 
Richter-Menge  (1991)  found  that  such  visual  esti¬ 
mates  maybe  off  by  ±5°.  This  error  was  found  after 
her  field  samples  were  sent  to  the  laboratory  for 
thin-section  analysis.  The  axial  DB  oc  values  vs. 
porosity  for  the  90°-0°  and  90°-90°  axial  DB  load 
tests  are  shown  in  Figure  17.  Also  shown  for  com¬ 
parison  is  the  curve  derived  from  the  equation 
proposed  by  Timco  and  Frederking  (1990)  for 
estimating  the  average  oc  strength  of  horizontally 


Table  4.  Test  data  for  first-year  sea  ice  loaded  perpen¬ 
dicular  to  the  c-axes  at  -16°C. 


Bulk 

Sample  density 
no.  (Mg/to5) 

Melt 

salinity 

(X.) 

Brine 

vol. 

(Xo) 

Air 

vol. 

(%c) 

Porosity 

(X.) 

Ice 

density  ac 
(Mg/m*)  (MPa) 

1 

0.864 

3.06 

11.3 

62.7 

73.9 

0.851 

5.12 

2 

0.889 

2.92 

11.1 

35.4 

46.5 

0.877 

6.47 

3 

0.885 

3.40 

12.8 

40.3 

53.2 

0.870 

6.40 

4 

0.897 

3.41 

13.0 

27.8 

40.8 

0.882 

5.73 

5 

0.892 

3.23 

12.3 

32.9 

45.1 

0.878 

6.83 

6 

0.894 

2.65 

10.1 

30.3 

40.4 

0.882 

5.72 

7 

0.891 

3.54 

13.4 

34.5 

48.0 

0.875 

6.51 

8 

0.890 

2.94 

11.2 

34.3 

45.5 

0.877 

6.49 

9 

0.864 

3.38 

12.4 

63.6 

76.0 

0.849 

4.62 

10 

0.896 

3.35 

12.8 

28.9 

41.7 

0.881 

6.11 

11 

0.887 

3.34 

12.6 

38.0 

50.6 

0.873 

6.10 

12 

0.896 

2.73 

10.4 

27.7 

38.1 

0.884 

6.55 

13 

0.885 

2.83 

10.7 

39.7 

50.3 

0.873 

6.57 

14 

0.888 

3.33 

12.6 

37.4 

50.0 

0.873 

5.57 

15 

0.894 

3.14 

12.0 

30.4 

42.4 

0.880 

6.44 

16 

0.899 

3.11 

11.9 

25.0 

36.9 

0.885 

5.90 

17 

0.886 

3.48 

13.1 

39.2 

52.4 

0.871 

5.12 

18 

0.885 

3.01 

11.4 

39.9 

51.2 

0.872 

5.49 

19 

0.884 

3.02 

11.4 

41.8 

53.1 

0.870 

6.55 

20 

0.8% 

3.12 

11.9 

29.0 

40.9 

0.882 

5.92 

21 

0.898 

2.92 

11.2 

25.9 

37.0 

0.885 

6.51 

Avg 

0.889 

3.14 

11.9 

36.4 

48.3 

0.875 

6.03 
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Table  5.  Test  data  for  first-year  sea  ice  loaded  parallel 
to  the  c-axes. 

Bulk  Melt  Brine  Air  Ice 

Sump/e  density  salinity  vol.  vol.  Porosity  density  ac 


no. 

(My/m3) 

<%o) 

<%°) 

(%.) 

(%.) 

(Mg/m3)  (MPa) 

-S°C  tests 

i 

0.908 

3.21 

31.4 

15.2 

46.7 

0.875 

5.78 

2 

0.893 

3.00 

28.8 

31.9 

60.8 

0.862 

6.43 

3 

0.893 

3.04 

29.2 

32.0 

61.2 

0.862 

6.06 

4 

0.900 

3.35 

32.5 

24.3 

56.8 

0.866 

6.04 

5 

0.903 

3.36 

32.7 

31.4 

64.1 

0.868 

4.49 

6 

0.907 

3.30 

32.2 

17.5 

49.7 

0.872 

6.80 

7 

0.908 

3.72 

36.4 

16.9 

53.2 

0.869 

6.14 

8 

0.903 

3.18 

30.9 

21.3 

52.2 

0.870 

6.02 

9 

0.901 

3.48 

33.8 

33.9 

67.7 

0.865 

4.53 

10 

0.903 

3.27 

31.8 

30.8 

62.6 

0.869 

4.36 

11 

0.900 

3.30 

32.0 

24.3 

56.3 

0.866 

5.58 

Avg 

0.902 

3.29 

32.0 

25.4 

57.4 

0.868 

5.66 

-16°C  tests 

1 

0.893 

3.50 

13.3 

32.2 

45.6 

0.877 

5.94 

2 

0.898 

3.70 

14.2 

26.7 

40.9 

0.882 

6.51 

3 

0.893 

3.63 

13.8 

32.7 

46.5 

0.877 

6.97 

4 

0.884 

3.33 

12.5 

42.1 

54.6 

0.869 

7.62 

5 

0.891 

3.18 

12.1 

33.7 

45.8 

0.877 

6.75 

6 

0.898 

3.74 

143 

27.3 

41.6 

0.881 

5.68 

7 

0.888 

3.46 

13.1 

37.6 

50.7 

0.873 

6.40 

8 

0.890 

3.41 

12.9 

35.5 

48.5 

0.875 

766 

9 

0.897 

3.63 

13.9 

27.5 

41.4 

0.881 

6.92 

10 

0.891 

3.39 

12.9 

33.9 

46.8 

0.876 

6.72 

11 

0.890 

3.61 

13.7 

35.8 

49.4 

0.874 

5.07 

12 

0.897 

3.08 

11.8 

27.1 

38.9 

0.883 

7.25 

Ayr 

0.892 

3.47 

13.2 

32.7 

45.9 

0.877 

6.63 

loaded  columnar  sea  ice.  Their  equation  for  this 
curve  is 

cxc  =  37(ef“[l-(M/270F5]  (3) 

where  n  is  porosity.  This  may  be  simplified  to 

oc  =  8.08  -  0.492(h)0  5  « 

for  a  strain  rate  e  of  10~3s_1.  Their  universal  equa¬ 
tion  relates  the  horizontally  loaded  uniaxial  un¬ 
confined  compressive  strength  explicitly  to  co¬ 
lumnar  sea  ice,  porosity  and  strain  rate,  and  im¬ 
plicitly  to  the  salinity,  temperature  and  brine-free 
density  of  the  sea  ice  (Timco  and  Frederking  1991). 
As  seen  in  these  figures,  the  axial  DB  ocvs.  n  trend 
is  similar  to  Timco  and  Frederking's,  but  the  axial 
DB  oc  values  are  higher.  This  may  be  due  to  the 
higher  apparent  strain  rate  under  which  our  tests 
were  made,  a  strain  rate  for  which  their  equation 
may  not  be  valid  (Timco  and  Frederking  1991). 
The  axial  DB  oc  results  may  also  be  compared  with 
the  first-year  sea  ice  uniaxial  oc  results  of  Richter- 
Menge  (1991),  which  were  made  at  -10°C  and  a 
strain  rate  of  10-3  sr1.  At  a  load  direction  of  90°-0° 


Table  6.  Test  data  for  first-year  sea  ice  with  aligned  c 
axes,  loaded  vertically. 

Bulk  Melt  Brine  Air  Ice 

Sample  density  salinity  vol.  vol.  Porosity  density  oc 


no. 

(Mg/m3) 

<%.) 

<%o) 

(Mg/m3) 

(MPa) 

i 

0.883 

3.25 

-S°C  tests 

30.9  43.3 

74.2 

0.850 

5.45 

2 

0.873 

3.10 

29.1 

53.3 

82.5 

0.842 

5.47 

3 

0.870 

3.07 

28.8 

56.5 

85.3 

0.839 

3.86 

4 

0.886 

3.65 

34.8 

40.4 

75.2 

0.849 

4.32 

5 

0.884 

3.32 

31.6 

42.0 

73.6 

0.850 

5.14 

6 

0.887 

3.45 

32.9 

39.2 

72.2 

0.851 

4.01 

7 

0.893 

3.73 

35.9 

32.8 

68.7 

0.855 

4.92 

8 

0.885 

3.94 

37.6 

41.4 

78.9 

0.845 

5.62 

9 

0.881 

2.81 

26.7 

44.5 

71.2 

0.852 

4.84 

10 

0.885 

3.04 

29.0 

40.5 

69.5 

0.854 

3.97 

11 

0.889 

4.16 

39.8 

38.2 

78.0 

0.846 

5.38 

12 

0.871 

3.09 

29.0 

55.3 

84.3 

0.840 

5.21 

13 

0.873 

3.10 

29.1 

54.0 

83.1 

0.841 

5.12 

14 

0.864 

3.76 

35.0 

64.5 

99.4 

0.826 

3.91 

15 

0.890 

2.97 

28.5 

34.7 

63.2 

0.860 

4.95 

16 

0.889 

3.44 

32.9 

36.6 

69.5 

0.854 

4.29 

Avg 

0.881 

3.37 

32.0 

44.8 

76.8 

0.847 

4.78 

1 

0.898 

3.11 

-7°C  tests 

22.4  26.1 

48.5 

0.874 

5.31 

2 

0.901 

2.87 

20.7 

22.9 

43.6 

0.878 

4.24 

3 

0.915 

2.88 

21.1 

6.8 

28.0 

0.892 

5.05 

4 

0.914 

3.75 

27.5 

9.9 

37.4 

0.884 

5.91 

5 

0.912 

3.91 

28.6 

11.5 

40.1 

0.881 

5.31 

6 

0.879 

4.29 

30.2 

47.7 

78.0 

0.846 

4.89 

7 

0.928 

3.79 

28.2 

5.7 

22.5 

0.897 

5.95 

8 

0.902 

3.66 

26.5 

22.1 

48.6 

0.873 

5.05 

9 

0.918 

3.99 

29.4 

5.2 

34.5 

0.996 

5.00 

Avg 

0907 

3.58 

26.1 

17.5 

42.4 

0.891 

5.19 

1 

0.920 

3.28 

-18°C  tests 

12.3  2.9 

15.2 

0.905 

5.25 

2 

0.908 

3.79 

14.1 

15.9 

30.0 

0.892 

5.79 

3 

0.908 

3.41 

12.6 

15.9 

28.5 

0.893 

6.28 

4 

0.922 

2.91 

11.0 

0.7 

11.7 

0.909 

7.42 

5 

0.917 

2.91 

10.9 

6.0 

16.9 

0.904 

7.14 

6 

0.907. 

3.97 

14.7 

17.6 

32.3 

0.890 

6.13 

7 

0.874 

4.69 

16.7 

53.8 

70.5 

0.855 

5.50 

8 

0.897 

3.14 

11.5 

27.7 

39.2 

0.883 

5.55 

9 

0.907 

4.06 

15.0 

18.1 

33.1 

0.889 

5.75 

Avg 

0.907 

3.57 

13.2 

17.62 

30.8 

0.891 

6.09 

she  reported  an  average  oc  of  7.30  MPa  for  sea  ice 
with  an  average  porosity  of  30.5  %o.  At  this  poros¬ 
ity  the  regression  curve  in  Figure  17a  gives  a  oc 
value  of  7.41  MPa.  For  her  tests  at  a  load  direction 
of  90°-90°,  the  average  porosity  and  oc  values 
were  33.0%o  and  6.56  MPa,  respectively.  This  oc 
value  compares  with  a  value  of  6.61  MPa,  which 
can  be  determined  for  the  same  porosity  from  the 
regression  curve  in  Figure  17b. 

The  regression  curves  were  used  to  determine 
that  the  oc  values  for  90°-0°  and  90°-90°  were  6.81 
and  6.32  MPa,  respectively,  at  n  -  40  %o  (Fig.  17). 
These  values  and  the  uniaxial  oc  values  of  Richter- 
Menge  (1991),  at  an  unknown  porosity,  are  plotted 


13 


3 


3 


a.  Loaded  parallel  to  the  predominant  horizontal 
c-axis  direction. 


b.  Loaded  perpendicular  to  the  predominant  hori¬ 
zontal  c-axis  direction. 


Figure  17.  Horizontal  first-year  sea  ice  uncon¬ 
fined  compressive  strength  vs.  porosity  at  a  strain 
rate  of  10~3  s_I . 


Figure  18.  Horizontal  first-year  sea  ice  un¬ 
confined  compressive  strength  vs.  the  angle 
between  the  applied  load  and  the  predomi¬ 
nant  c-axis  direction. 
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in  Figure  18.  Also  shown  is  the  horizontal  oc  trend 
vs.  the  angle  between  the  applied  load  and  the  pre¬ 
ferred  c-axis  direction  as  estimated  by  Wang  (1979) 
from  tests  made  at  -10°C  and  a  strain  rate  of  10~3 
s'1  and  the  trend  established  by  Peyton  (1966) 
from  tests  at  a  slower  strain  rate,  which  he  gives  in 
terms  of  a  load  rate  of  -0.053  MPa/s.  The  temper¬ 
ature  at  which  his  tests  were  made  was  not  given 
in  his  report. 

The  horizontal  uniaxial  oc  anisotropy  trend 
shown  in  Figure  18  is  related  to  the  sea  ice  struc¬ 
ture  and  is  similar  to  the  strength  anisotropy  trends 
observed  in  the  data  from  uniaxial  oc  tests  made  on 
bedded  rock  (Brown  et  al.  1977).  The  trend  estab¬ 
lished  by  Peyton  (1966)  was  well  determined 
through  a  large  number  of  uniaxial  <jc  tests.  The 
curve  proposed  by  Wang  (1979)  was  less  well 
determined  experimentally.  However,  it  can  be 
inferred  from  the  tests  of  Richter-Menge  (1991),  as 
plotted  in  Figure  18,  and  the  data  shown  in  Borod¬ 
kin  et  al.  (1992)  that  further  testing  is  required  to 
clarify  the  strength  anisotropy  trend  indicated  by 
the  Wang  and  Peyton  curves  in  Figure  18.  The 
axial  DB  oc  values  at  ti  -  40  %o,  shown  in  Figure  18, 
are  comparable  with  the  uniaxial  ac  values  of 
Richter-Menge  but  do  not  reveal  the  magnitude  of 
the  strength  anisotropy  suggested  by  the  Wang 
curve.  Sinha  (1983)  also  observed  little,  if  any, 
difference  between  the  uniaxial  oc  values  from 
tests  made  parallel  vs.  perpendicular  to  the  pre¬ 
ferred  c-axis  direction.  His  tests  were  made  at 
-10°C  and  a  strain  rate  of  about  -10-4  s'1. 

Axial  DB  load  tests  were  made  on  Stefanson 
Sound  first-year  sea  ice  with  the  load  applied  in 
the  growth  direction,  parallel  to  the  columnar 
crystal  axes.  Tests  were  made  at  -5°,  -7°  and  -18°C 
as  listed  in  Table  6.  The  average  oc  values  at  -7° 
and  -18°C  were  5.19  and  6.09  MPa,  respectively. 
The  axial  DB  <jc  difference  represents  a  tempera¬ 
ture  effect  of  0.082  MPa/°C. 

Another  set  of  vertical  axial  DB  tests  were  made 
at -7°  and-17°C  on  first-year  sea  ice  obtained  from 
a  melt  pool  in  a  multiyear  floe  (Table  7).  In  this  sea 
ice  the  crystals  had  an  average  width  of  10  mm, 
and  the  c-axes  did  not  have  a  preferred  horizontal 
alignment.  The  difference  in  strength  at  each  tem¬ 
perature  indicates  a  oc  variation  of  0.07  MPa/°C. 
Note  that  the  porosity  at  -7° C  is  about  66%  more 
than  the  value  at  -17°C.  This  is  related  to  the 
change  in  brine  volume  that  occurs  with  tempera¬ 
ture.  When  the  temperature  of  the  sea  ice  decreas¬ 
es,  a  portion  of  the  water  in  the  brine  inclusions 
freezes  to  the  chamber  walls.  This  occurs  in  order 
to  maintain  the  brine  concentration  in  thermal 


equilibrium.  When  the  ice  warms,  the  reverse 
occurs;  that  is,  brine  dilution  occurs  as  the  en¬ 
riched  brine  causes  the  chamber's  freshwater  ice 
walls  to  melt. 

At  -7°  and  -18°C  the  average  oc  value  for  the  c- 
axis  aligned  Stefanson  Sound  sea  ice  was  5.19  and 
6.09  MPa,  respectively.  For  the  non-aligned  melt 
pond  ice  the  average  oc  value  was  5.21  and  5.91 
MPa  at  -7°  and  -17°C,  respectively.  Because  these 
values  are  in  good  agreement,  it  would  appear 
that  the  vertical  ac  strength  of  columnar  sea  ice  is 
not  affected  by  horizontal  c-axis  alignment. 

All  the  vertical  axial  DB  oc  values  are  plotted  vs. 
n  in  Figure  19.  At  a  porosity  of  40 %o  the  regression 
curve  through  the  data  gives  a  oc  of  5.44  MPa.  This 
value  is  low  when  compared  to  the  horizontal 
axial  DB  load  test  results  in  Figure  17,  and  it  is  not 
in  conformity  with  vertical  vs.  horizontal  uniaxial 
unconfined  compression  tests  made  on  first-year 
sea  ice. 

To  illustrate  the  latter,  Sinha  (1983)  reported 
vertical  uniaxial  ac  values  about  2-5  times  higher 

Table  7.  Test  data  for  melt  pond  sea  ice,  with  random 
crystal  c-axis  orientations,  loaded  vertically. 

Bulk  Melt  Brine  Air  Ice 


Sample 

no. 

density 

(Mg/m3) 

salinity 

<%o) 

i  ml. 

(%o) 

vol. 

(%,) 

Porosity 

(%°) 

density 

(Mg/m3) 

Oc 

(MPa) 

-7°C  tests 

i 

0.910 

4.31 

31.4 

14.4 

45.8 

0.876 

5.33 

2 

0.908 

4.50 

32.7 

17.2 

49.9 

0.872 

5.22 

3 

0.912 

3.98 

29.1 

12.0 

41.1 

0.880 

5.53 

4 

0.910 

4.88 

35.6 

15.7 

51.2 

0.871 

5.92 

5 

0.911 

4.00 

29.2 

13.2 

42.4 

0.879 

5.51 

6 

0.909 

4.12 

30.0 

15.3 

45.3 

0.876 

4.26 

7 

0.907 

4.58 

33.3 

17.8 

51.1 

0.871 

4.79 

8 

0.909 

4.27 

31.1 

15.9 

47.1 

0.875 

5.48 

9 

0.907 

4.58 

33.3 

18.1 

51.4 

0.871 

5.18 

10 

0.894 

3.68 

26.4 

31.4 

57.7 

0.865 

5.53 

11 

0.909 

4.57 

33.3 

16.1 

49.4 

0.873 

4.74 

12 

0.911 

4.02 

29.4 

13.2 

42.6 

0.879 

4.94 

13 

0.910 

3.50 

25.5 

13.7 

39.2 

0.882 

5.33 

Avg 

0.908 

4.23 

30.8 

16.5 

47.2 

0.875 

5.21 

-17°C  tests 

1 

0.910 

4.19 

15.6 

15.1 

30-6 

0.891 

5.47 

2 

0.980 

3.45 

12.8 

15.6 

28.4 

0.893 

5.28 

3 

0.914 

4.04 

15.1 

9.7 

24.8 

0.897 

7.65 

4 

0.899 

3.65 

13.4 

26.0 

39.4 

0.883 

6.95 

5 

0.915 

5.09 

19.0 

9.7 

28.8 

0.893 

5.64 

6 

0.909 

3.51 

13.0 

14.5 

27.6 

0.894 

5.34 

7 

0.908 

3.62 

13.4 

15.7 

29.1 

0.893 

5.78 

8 

0.900 

3.41 

12.5 

24.3 

36.8 

0.886 

5.58 

9 

0.902 

3.89 

14.3 

22.9 

37.2 

0.885 

5.95 

10 

0.907 

3.33 

12.3 

16.6 

28.9 

0.893 

5.12 

11 

0.911 

3.99 

14.9 

13.0 

27.9 

0.894 

6.17 

12 

0.907 

4.27 

15.8 

18.1 

34.0 

0.888 

5.86 

13 

0.910 

4.31 

16.0 

14.4 

30.4 

0.891 

6.01 

Avg 

0.913 

3.90 

14.5 

16.6 

31.1 

0.891 

5.91 
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n,  Porosity  (%.)  \  q-3  S-1 

than  those  obtained  from  horizontal  tests  at  -10°C  plane,  parallel  to  the  crystal  lattice  basal  plane  and 

and  a  strain  rate  between  about  4  x  10"5  and  7  x  the  columnar  ice  structure.  Therefore,  it  is  most 

1CT5  s'1.  Frederking  and  Timco  (1989)  observed  a  reasonable  that  the  axial  DB  test  would  give  the 

uniaxial crc  value  three  times  larger  in  their  vertical  results  noted. 

vs.  horizontal  load  tests  at  -1 1°C  and  at  strain  rates 

between  2  x  IQ"5  and  5  x  10-4  s'1.  Kuehn  and  Schul-  DISCUSSION 

son  (1993)  made  uniaxial  tests  on  laboratory-grown 

columnar  saline  ice  with  unaligned  c-axes  in  the  The  axial  DB  test  results  presented  in  this  report 

horizontal  plane.  The  only  consistency  shown  over  indicate  that  this  test  is  well  suited  for  determining 

a  temperature  range  of -5°  to  -40°C  and  strain  rates  the  unconfined  compressive  strength  of  multiyear 

from  10*6  to  10"1  s'1  was  that  the  vertically  loaded  sea  ice  and  freshwater  ice.  The  agreement  between 

ice  failed  at  a  stress  1. 1-7.5  times  higher  than  the  the  axial  DB  and  laboratory  uniaxial  unconfined 

horizontally  loaded  samples.  At  a  strain  rate  of  compression  tests  made  on  these  ice  types  was 

10"3  s'1  and  temperatures  of  -5°,  -10°,  -20°  and  found  to  be  excellent  at  a  strain  rate  of  10"3  s'1. 

-40°C,  the  vertically  loaded  samples  were  1.2, 2.8,  The  axial  DB  tests  on  first-year  sea  ice,  a  less 
1.1  and  2  times  stronger,  respectively,  than  the  "brittle"  material,  were  more  difficult  to  assess 

horizontally  loaded  samples.  Richter-Menge  (1986)  because  of  the  inherent  variation  in  sea  ice  struc- 

showed  first-year  sea  ice  to  be  about  1.5  times  hire  and  physical  properties.  Numerous  strength 

stronger  when  loaded  vertically  to  the  crystal  col-  tests  on  sea  ice  need  to  be  made  to  establish  trends 

umns  than  when  loaded  horizontally.  Her  uncon-  related  to  crystal  size,  load  orientation  and  rate  vs. 

fined  compression  tests  were  made  at  90°  to  the  the  ice  structure,  ice  porosity,  temperature,  etc. 

mean  c-axis  direction  at  a  test  temperature  of  -1 0°C  The  axial  DB  horizontal  load  test  results  did  show 

and  a  strain  rate  of  lO""3  s"1.  However,  Schwarz  a  oc  dependence  on  the  direction  of  the  applied 

(1970)  found  that  the  uniaxial  oc  of  Baltic  Sea  ice,  load  relative  to  the  c-axis  alignment.  The  results 

tested  at0°,-10°  and-20°C  andatstrain  rates  from  also  appear  to  give  a  good  assessment  of  the  mean 

3  x  10"4  to  3  x  10"1  s'1,  was  always  stronger  when  horizontal  sea  ice  strength,  which  is  the  preferred 
loaded  in  the  horizontal  vs.  vertical  direction!  strength  for  use  in  the  design  of  offshore  structures 

Except  for  the  contradictory  results  of  Schwarz,  (W ang  1979). 

there  is  a  strong  indication  that  sea  ice  fails  at  a  In  this  report  a  <TC  temperature  correction  factor 

higher  uniaxial  oc  when  the  ice  is  loaded  in  the  of  0.075  MPa/°C  was  used  to  correct  freshwater 

growth  direction  vs.  parallel  to  it.  However,  it  is  not  and  multiyear  sea  ice  test  results  to  a  common 

possible  at  this  time  to  quantify  what  the  vertical-  temperature.  This  value  was  reported  to  be  in 
to-horizontal  strength  ratio  is  for  a  given  set  of  sea  agreement  with  the  results  of  other  authors  and 
ice  properties  and  loading  conditions.  With  regard  was  found  to  be  in  agreement  with  the  test  results 

to  the  anomalous  axial  DB  vertical-vs.-horizontal  reported  here.  However,  this  temperature  correc- 

load  test  results,  it  can  be  expected  that  ball  loading  tion  is  not  in  agreement  with  the  test  results  of 

parallel  to  the  crystal  columns  will  induce  stresses  Carter  (1970),  Schwarz  (1970),  Haynes  (1979)  and 

that  result  in  the  ice  splitting  along  its  weakest  Schulson  (1990).  These  investigators  found  tem- 
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perature  correction  factors  for  uniaxial  unconfined 
compressive  strength  to  be  two  to  seven  times 
higher  than  the  one  used  in  this  report.  Further 
work  is  needed  to  resolve  this  inconsistency,  as  well 
as  the  issue  of  sea  ice  strength  vs.  porosity  and 
temperature.  For  example,  a  given  sea  ice  porosity 
can  be  found  at  different  temperatures.  Therefore, 
all  other  conditions  being  equal,  the  ice  strength 
will  be  temperature  dependent.  The  strength  equa¬ 
tion  of  Timco  and  Frederking  (1990)  presented  in 
this  report  does  not  account  for  this  porosity-vs.- 
temperature  effect,  which  must  be  responsible  for 
some  of  the  scatter  observed  in  the  unconfined 
compression  test  results. 

Another  issue  not  addressed  in  this  report  is  the 
effect  of  sample  geometry  and  size  and  the  various 
types  of  end  conditions  that  may  have  affected  the 
uniaxial  oc  values  obtained  by  other  investigators. 
For  example,  "bonded  end  caps  exert  a  lateral  con¬ 
finement  which  suppresses  both  crack  nucleation 
and  the  propagation  through  the  end  zones  of  axial 
splits  which  initiate  during  loading"  (Schulson 
1990).  This  confinement  results  in  failure  strengths 
on  the  order  of  20%  higher  than  if  the  ends  were  not 
constrained.  Such  test  confinement  vs.  non-con¬ 
finement,  as  well  as  sample-machine  alignment 
and  different  sample  slenderness  ratios,  may  have 
contributed  to  the  variation  in  the  uniaxial  ac  val¬ 
ues  reported  in  the  literature  and  used  for  compar¬ 
ative  purposes  in  this  report. 

The  axial  DB  test  will  not  provide  ac  values  at  a 
strain  rate  other  than  10~3  s-1,  nor  will  the  test  pro¬ 
vide  stress-strain  information  for  determining,  for 
example.  Young's  modulus.  However,  a  simple 
sonic  velocity  meter,  as  used  in  this  study,  may  well 
serve  this  purpose.  In  rock  mechanics  it  has  been 
shown  that  Young's  modulus  can  be  well  correlat¬ 
ed  with  the  dynamic  modulus  derived  from  sonic 
velocity  tests  (Eissa  and  Kazi  1988).  With  additional 
testing  of  ice  a  similar  correlation  could  be  found. 

In  rock  mechanics  the  point-load  strength  fre¬ 
quently  replaces  the  uniaxial  unconfined  compres¬ 
sive  strength  because  the  test  is  quicker  to  perform 
and  less-stringent  sample  preparations  are  required . 
The  test  may  be  readily  done  in  the  field  at  nearly 
in-situ  conditions  and  has  been  shown  to  be  as  reli¬ 
able  as  the  uniaxial  unconfined  compression  test, 
with  its  expensive  electrohydraulic  loading  and 
sample  preparation  systems.  Similar  exacting  sam¬ 
ple  preparation  procedures  and  testing  equipment 
are  required  in  the  uniaxial  unconfined  compres¬ 
sion  testing  of  ice.  For  axial  DB  tests  the  only  equip¬ 
ment  required  is  a  simple  loading  jig,  a  hand- 
operated  hydraulic  pump  ram  and  a  load  cell  with 


a  digital  peak-load  indicator.  Auxiliary  equipment 
would  include  a  cutoff  saw,  used  to  section  core 
into  lengths  of  1 .06  ±  0.01  times  the  core  diameter, 
and  appropriate  measuring  equipment  for  deter¬ 
mining  sample  density  and  temperature.  The  find¬ 
ings  of  this  report  indicate  that  the  axial  DB  test  is 
ideally  suited  for  determining  the  unconfined  com¬ 
pressive  strength  of  ice  in  the  field. 
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